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Abstract: We present a new observable sensitive to a certain combination of CP violating 
phases in supersymmetric extensions of the Standard Model, viz. a triple product of momenta 
in the cascade decay of a heavy squark via an on-shell neutralino and off-shell slepton. We 
investigate the regions of parameter space in which the signal is strong enough to be detectable 
at the LHC with ~ (lO'^ — 10^) / sin^(2A(/j) identified events, where Aip is a certain combination 
of phases in the MSSM presented in the text. 
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1. Introduction 



In the past decade it became clear that the Standard Model cannot account for electroweak 
baryogenesis (for recent reviews, see 0, The CP violating phases (the CKM phase and 

the strong phase) are too small, and with the Higgs mass constrained by LEP2 to be mn > 
115 GeV the electroweak phase transition is not strong enough to suppress the sphaleron 
process. 

However, it is still possible that baryogenesis is connected with electroweak scale physics 
and has to do with the other missing ingredient from the Standard Model, which is the mech- 
anism for stabilizing the electroweak scale. If that mechanism is supersymmetry, then many 
new possibilities open up. The new supersymmetric sectors contain, in principle, many new 
phases. Also, there are many more contributions to the Higgs potential, and the electroweak 
phase transition can be made stronger. 

In the context of the MSSM, there is only a small window in parameter space left to 
accommodate electroweak baryogenesis §, |^ §• One difficulty arises because there are 
no cubic terms in the tree-level Higgs potential, and the bounds from LEP2 on the Higgs 
and stop masses severely restrict the loop-induced contributions and disfavor a strong first- 
order phase transition. Nonetheless, there seems to be no difficulty in generating a strong first 
order transition in slightly more general supersymmetric models, such as those involving extra 



singlet Higgs in which there can be cubic terms at tree-level y, 10, 11 , |l^, 13, 14 , |15| . The other 
difficulty is in ensuring sufficiently large CP violating phases, given the stringent bounds on 
certain combinations of such phases from electron and neutron electric dipole moment (EDM) 
experiments [|l6|, 0, |l|, |l9|. Since the operator responsible for EDM involves a helicity flip, 
these bounds mostly constrain the phase combination involving the //-parameter and the 
wino M2 phases. The same phase combination is responsible for the dominant electroweak 
baryogenesis mechanisms in the MSSM, leading to tension with the EDM constraints. This 
may also be relaxed in extended models involving additional phases [p|-[|T^],[20|. 

Motivated by the possible abundance of CP violating phases in supersymmetric extensions 
of the Standard Model ||2l|], their relevance to baryogenesis, and their impact on Higgs and 
sparticle spectra p2| , p^ j it is important to explore them in a variety of ways. In particular, 
there have been a number of suggestions for their direct detection in collider experiments 



|2g, g6|, 123, g|, |3g, |3g]. 

The MSSM soft lagrangian contains many different phases, however, not all are physical. 
As emphasized in |3^, the physical phases are those which are invariant under both Uji{l) 
and Upq{1) symmetries. There are several obvious such invariant: (i) the phases in the 
off diagonal elements of the soft scalar masses; (ii) the relative phases between the gaugino 
masses; (iii) the relative phases between the different A parameters. The phases which are 
affected by the R and PQ rotations are 0^, 0f, (the phases of the fi and b terms), (pMa 
and (t>^^, the overall phases of the A parameters. One can then build reparameterization 
invariant linear combinations of these phases which appear in physical processes (for one 



such parameterization see p4[| ) 
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EDM experiments are sensitive to linear combinations of phases involving (p^ since the 
dominant diagram involves the chargino exchange (see for example []T8| ). In this paper we 
present a new observable which is sensitive to a different combination of the phases than is 
probed in the EDM experiments. It involves the phases in the couplings of the neutralinos 
to the sleptons, and in a typical limiting case depends on the difference between the phases 
of the bino and wino mass parameters Mi and M2 (which can differ for nonuniversal gaugino 
masses). Since it requires no higgsino insertion it is insensitive to the phase of the //-parameter. 

We propose an asymmetry parameter related to the usual, T-violating, triple product 
{pi ■ (p2 X ps)), with Pi being 3 independent vectors in a given reaction. This quantity 
constitutes a direct measurement of T-invariance violation, which, in the context of CPT 
invariant theories, is a measurement of CP-violation. 

As can be expected on general grounds, any such asymmetry parameter is the result 
of interfering diagrams. As such, it is usually suppressed with respect to the leading non- 
interfering contribution. In the case of a reaction which proceeds through an on-shell cascade 
decay, interference terms are suppressed with respect to the on-shell amplitudes by a factor of 
the width. However, if some part of the reaction is forced to proceed off-shell, the interference 
terms are comparable to the leading amplitudes. We exploit this fact in the decay of a squark 
into a quark and two leptons q ^ q + N2 q + l'^ +l~ + Ni, where are neutralinos and 
we assume that A^i is the lightest supersymmetric particle (LSP). We show that in regions 
of parameter space where the decay of is through an off-shell slepton the asymmetry 
parameter, r/ (defined below), can be very large, rj ~ 0{1). 

The number of events required for determining r/ scales as oc l/rf (assuming Gaussian 
statistics). This calls for as precise a theoretical estimate of r/ as possible, since a factor of 3 
reduction can translate into an order of magnitude more events. Clearly, the actual number 
of events needed is affected by many experimental considerations as well. Reduction of the 
signal due to mis-tagging, detector resolution, tagging efficiency, etc., will inevitably increase 
the number of events required. After presenting the theoretical results we attempt to estimate 
the signal reduction due to experimental limitations and present the number of events needed 
to determine the existence of a CP asymmetry. 

The paper is organized as follows: In section ^ we present a general discussion of the T, 
CP, and P-violating triple product and the relevant observables. In section |3| we consider the 
reaction t ^ t + + l~ + Ni and show that it contains appropriate observables sensitive 
to certain CP-phases in the MSSM. However, we find that if the reaction proceeds through 
an on-shell cascade decay, the signal is too small to be measured. It is also possible for the 
reaction to proceed via an off-shell cascade. In that case, as we shall see in section ^, the 
effect is greatly enhanced and the signal may be large enough to be detectable. In section ^ 
we attempt to estimate the reduction in the signal strength due to experimental limitations 
and section |6| contains our conclusions. 
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2. General discussion 



In this paper we will concentrate on CP-violating effects present in cascade decays of heavy 
supersymmetric particles. When an unstable particle decays through a reaction involving 
4 other particles, it is possible that the expression for the rate contains a contribution of 
the form^ ^^i/o/JPoPiPf^'s ' where pi are four independent momenta. In the rest frame of the 
decaying particle, po = (Mq, 0,0,0), this term is the usual P and T-odd observable given by 
the expectation value of the triple product 

r = -Mo Pi • {p2 X ps) , (2.1) 



where Mq is the mass of the decaying particle. As is well known ||35[], any measurement of a 
non-zero expectation value (T) implies both P violation and either T violation or a "strong 
phase". Assuming that CPT invariance is unbroken, T violation is equivalent to the non- 
conservation of CP, and is manifested by a CP-violating phase in the Lagrangian. A strong 
phase refers to a CP-conserving phase, due, e.g., to a strong or electromagnetic final state 
interaction or the phase associated with the width in the propagator of an unstable particle. 
The two effects can in principle be separated by measuring the expectation values of both T 
and T, where T = —M^pf ■ {pf x pf) is the corresponding triple product for the decay of 
the antiparticle, and pi are the physical momenta of the antiparticles in the final state. Since 
T is both P and T-odd a non-zero (T) requires the interference of two contributions to the 
amplitude involving different phases and different parity, i.e., the amplitude in the rest frame 
Po = must contain 

{pi\H\pQ) = A{p{)e'^P^+'f'^'^ + B{pi)e'^^+P'^+'t"^\ (2.2) 

where pi refers collectively to the final momenta, and, 

A{-p) = A{p) (2.3) 
B{-p) = -B{p). 

(j)A,B and pA,B are respectively the so-called weak and strong phases, which do (do not) 
change sign in the CP-conjugate process^. The phase 7r/2 in the second term could have been 
absorbed m ps, but is instead pulled out for convenience. It always occurs in the relevant 
interference term between the parity even and odd amplitudes for (T) when one sums over 



the spins [|35[. For the example considered in this paper, it is just the explicit factor of i 
occurring in the trace of ^^'^^"jvlalp- Using ( ^.2]) , it is then straightforward to show that 



(T) = 2/C [sin(pA - pb) cos{<I)a - (Pb) + cos(pA - Pb) sin(0A - (Pb)] , (2.4) 

where /C is proportional to the phase space integral J dPS{T AB). The first (second) term 
in ( |2.4| ) requires a nonzero difference between the strong (weak) phases. Using CP, the 



^Our conventions are g^v = diag(l, —1, —1, —1) and eoi23 = — 1- 

^For a more thorough discussion of these phases see the excellent textbooks |3^, 
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corresponding amplitude for the antiparticle decay is 

{p^\H\pQ^) = A{-pt)e^'^P^-'t>^^ + B{-pt)e'^^+P^-'t>^\ (2.5) 



up to an irrelevant overall phase associated with our CP conventions. Comparing with ( |2.2| ) 
and using the parities of A and B in ( |2.3D , one sees that (T) differs from (T) by an overall 
sign (due to the fact that the observable is P-odd) and by (j)A,B —4>A,B, 

(T) = -2/C [sin(pA - Pb) cos((/)a - 4>b) - cos{pA - Rb) sin(0yi - ^b)] . (2.6) 

In particular, the T-odd term can be isolated by summing the particle and antiparticle asym- 
metries 

(T) + (T) = 4/C cos(/>A - pb) sin(</>A - (t>B). (2.7) 

In this paper we will present a manifestation of these general considerations in the par- 
ticular super symmetric cascade decays 

i^t + Na^t + l+ + 1- + Ni, (2.8) 

and their CP conjugates ^ + Na + l~ + + Ni, where Na (which we usually take 
to be A'^2) is assumed to be on-shell. These reactions have enough independent momenta and, 
as we shall show below, lead to non-vanishing expectation values for the triple product, Eq. 
( |2.lD , made of the top and di-lepton momenta. T and T refer respectively to the observables 
-M^ pt ■ {pi+ X pi-) and -M^ ptc ■ {pi- x pi+). 

It proves useful to formulate the discussion in terms of a dimensionless parameter em- 
bodying the CP asymmetry. This parameter is closely related to the triple product with the 
added advantage of allowing for a straightforward evaluation of the number of events needed. 
As shown in Fig. ||, in the rest frame of the on-shell Na, the incoming i and outgoing top 
define a z-axis. Momentum conservation forces the outgoing anti-lepton, lepton and the LSP 
to define a plane. A non-zero expectation value of T implies a non-zero average angle between 
the plane and the z-axis. We therefore define the asymmetry parameter 

_N+-N^ _N+-N_ 

where 

/■I dV 1'^ dV 

N+ = / -dcos^, iV_ = / -dcosO. (2.10) 

Jo a cos a J_iacosU 

Below, we derive an exact expression for r] in the neutralino's rest frame. Clearly, r] 
is not a relativistically invariant variable. If the LSP escapes detection it is impossible to 
reconstruct the neutralino's rest frame. Therefore, r] can only be constructed in the lab frame 
(i.e., detector frame), which inevitably affects the signal. In section ^ we show that this lack 
of knowledge of the correct frame gives rise to a dilution factor D, similar to that encountered 
in B physics experiments, 

Vexp = 'Dtith- (2.11) 
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Figure 1: The reaction geometry in tlie rest frame of Na ioi t t + + I + Ni. 9 is the angle 
between pt and pi+ x pi- . For — > + + Z+ + TVi, is the angle between pt^ and pi- x pi+ . 

Assuming Gaussian statistics, the number of events needed to make a statistically signif- 
icant measurement is given by 

^ = J- = (2-12) 

'lexp ^ 'Ith 

The dilution factor, as its name implies, leads to an increase in the number of events needed. 
In what follows we evaluate the asymmetry parameter "qth and show that in certain kinematical 
regimes it can be rather large 0(1). Also, we find that the dilution factor T) need not be very 
small and in fact does not present a very serious obstacle. 



3. CP-violation in the stop cascade decay via an on-shell slepton 

In this section we compute the expectation value of ^^uapPiPtP'i+P^- fo^' the cascade decay 
i ^ t + 1+ + I- + Ni via the two diagrams shown in Fig. ^ and written explicitly in the 
appendix. 

The relativistically invariant expectation value is given by 

/ M . a /3 \ Pj PtP?+ pf- . 

Pi PtPi+ Pi- > = rjr ' (^-^^ 



JdV 
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Figure 2: The Feynman diagrams corresponding to the matrix elements iMa (left) and iMb (right). 
The full matrix element is iM = iAia + iMb- 

where p^, pt, Pi+ and pi- are the four-momenta of the squark, quark, anti-lepton and lepton, 
respectively, and dT is the differential decay width. 

Assuming that the neutralino Na is on-shell the expectation value can be evaluated ex- 
actly; the details are presented in the appendix. Considering only flavor diagonal interactions, 
the numerator in Eq. ( ^.1| ) is 

dVe^^ap p- p, p,^ p^_ - 3 256^3 (^7— ) [J ^ ) ^^'^^ 

x{\9ir-Kr) 



Ja* \ / „(1 



9T 



where 



\Pt\' = ^ {{M? - Ml - Mff - AMlM^) . (3.3) 
The dimensionless function /(x+,x_) is given by 



/(j;+, x_) = {1 — Hi — x+ — X- + x^x-) {x+ + X- + fii — 1) (3-4) 



[1 - 111 - x+Y + //?7;-) {{I - /i; - x_)2 + ^2^;- 



where 

Mf Ml r? 

Na Na I 
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To arrive at an expression for the expectation value we need the total width of the 
reaction. In the narrow-width limit the computation is straightforward and can be carried 
out analytically. The result is given by 

^{Hi-tiif (3.6) 
n 




{Ml - Ml - Ml) {\gir + Ig'^r) + ^MtM^ReisTgr)) • 

Clearly, the general expression for the expectation value is fairly complicated. It is 
instructive to consider some limits where it simplifies. For example, if we assume there is no 
suppression from the kinematical factors (e.g. (1 — iJ,j) ~ 1) in Eq. (^.6|), and considering the 
limit where Na is a pure wino and A'^i a pure bino, the final expression is, 

(W Pt P?+ Pf- > = VJ^i J dx+dx.fix+,x.) (3.7) 

X sin (2ip'l - 2ip'l) 

where we have written the complex couplings as (7^ = \g^l \ exp(i(^^), etc. This expression is 
suppressed by the width of the slepton and one additional phase-space factor. This is simply 
a reflection of the fact that the signal is the ratio of an off-shell process to an on-shell one, 
i.e., the two diagrams in Fig. |2| cannot simultaneously be on-shell except for a set of measure 
zero. The integral over /(x+, x_) is of order unity and cannot enhance the signal. If the LSP 
is mostly a bino then the slepton decay width is roughly 

1 F;- 1 

~ — ~ . (3.8) 

7rM|- TT 300 ^ ' 



In order to reliably estimate the number of events needed to reach experimental sensitivity 
one must form a dimensionless quantity, such as the asymmetry variable presented in the 
previous section, Eq. (2.9). The dimensionful phase-space factors in Eq.(|3.7|) roughly cancel 



out in such an observable. Therefore we expect that the quantity in (^) gives us a good order 
of magnitude estimate for the number of events needed. Without even taking experimental 
limitations into account we need at least 10^ events to reach statistical significance. 

However, if one is close to the decay threshold such that M?^ — < ^i^v decay 
rate is suppressed and the asymmetry is enhanced. A more likely possibility is a spectrum 
where the slepton is forced to be off-shell. In this case there is no width suppression. We 
explore this possibility in the next section and show that indeed the signal is greatly enhanced. 
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4. Stop cascade decay via an off-shell slepton 



M^^ > Mz then the Z is on-sheh and 
M^^ < Mz then the Z 



In this section we consider the case where > . The neutrahno may decay through the 
3-body channel Na l'^ + + Ni via an off-shell slepton. While most of the results hold 
for general Na-, to simplify the discussion we will often take Na to be approximately wino. 
In general, there may be additional decay paths to consider, and in particular the neutralino 
can decay directly into A'^i and a Z-boson. If 
this channel dominates over the 3-body mode. However, if 
is off-shell and this reaction might compete with the diagram involving an off-shell slepton. 
Which is dominant is a detailed question depending on the spectrum. The coupling NaNiZ 
is a result of mixing with the higgsino component. For a = 2 it is therefore governed by the 
size of the n term compared with the gaugino masses Mi and M2 and is decoupled in the 
limit of large n. Our major goal is to illustrate the possibility of measuring a CP-violating 
effect rather than to exhaustively examine all of parameter space. In what follows we will 
therefore ignore the possible contribution of the NaNiZ vertex. In the more general case the 
interference with the Z diagram could enhance or reduce the effect. 

We expect the signal to have no parametric suppression as in the case of an on-shell 
decay discussed in the previous section. Also, to evaluate the number of events needed we 
concentrate on the asymmetry parameter r] defined in Eq. ( p. 91 ). The details are very similar 
to the previous section except that the interference terms in the width cannot be neglected. 
Therefore, for integrated luminosity C the total number of events is given by 



c 



dPSj 

2Mr 



\Ma? + \Mb\'' + 2Re{MaMl)) . 



(4.1) 



The evaluation of the phase-space integrals is presented in the appendix. The difference 
between the number of events in the upper and lower hemispheres is given by 



N. 



c 




\9l?) 



2Ml lTa(g'^*g'r9R9L 



I 



X / dxj^dx^ 



{1 — — Xj^ — X- + Xj^X_) {xj^ + X- + — 1) 



1/2 



(1 - x+ - (1 



(4.2) 
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The integrals evaluate to 

x[{m?- Ml - M^) {\gir + + ^MtM^^ Re {gfgT 

'l-Mi ^2n-X-/Ui)2 1 

X / dx- 







1 — X (1 — X — //;-)^ 



and 

dPSi 




2M, = I -F^ I ( / ^ ) (4-4) 



M| - Ml - M?) {\gir + Kl') - ^M,M^^ Re {g^gr 
dxj^dx^ - 

(/zi + x+ + x_ - 1) Re [g'lg'lgtat + sVhtdR) 
- 2(1 - x+ - x_ + - /ii)Re (ffi^5L55?*55?*)) • 

The entire expression is relativistically invariant, except for the limits in Eq. ( p.l[l| ) used 
to derive Eq. ( f4.2D , which are computed in the rest frame of Na- However, since (-^-) 
involves an integration over the entire upper (lower) hemisphere, these expressions are still 
invariant under boosts in the stop's direction which do not flip the direction of the top. In 
particular, the asymmetry parameter is unmodified when boosting to the rest frame of the 
stop. This is an important fact. It implies that the signal one constructs in the lab is only 
degraded by one's ignorance of the initial boost of the stop in the lab frame. 

In the case of an off-shell slepton the asymmetry variable r/ is an 0(1) number. The exact 
expression is given by the ratio of Eq. (O) to Eq. (O). There are several limiting cases 



where the final result is extremely simple. In particular, in the case where Na is a pure wino 
and A^i is a pure bino, assuming there are no strong kinematical suppressions and /i[ ^ 1, 
the expression simplifies to 

2 VGi(/ii) + G2(/ii)cos(2A99)y 

where we expressed the complex couplings as = exp{iip^£) and 

= - ^^£. (4.6) 

In the approximation of ignoring slepton mixings this is just the difference between the original 
phases of the neutralino masses before they were absorbed into the couplings. The kinematic 
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functions in Eq. (^]^) are given by 



G'2(/il] 



1 — /ii — x+ — x_ + x+x„) (x+ + X- + — 1] 



1/2 



(1 - /^i)(l + lO/ii + /i?) + 6/ii(l + /ii) log(/xi) 



In Fig. |3|we plot the asymmetry parameter rj vs. /ii for several choice of As claimed 
above, when the slepton is off-shell the asymmetry r/ can be very large, proportional to the 
phase times an 0(1) number as shown in Eq. ( f4.5[ ). For example, if we take Mfj^/M^ > 0.7 
we find that, ignoring experimental limitations, the number of events needed to make a 
determination of CP-violation in this cascade decay is approximately 



100 



sin^(2Av?) 



(4.7) 



In Table |T| we present the iit'j^ production cross-section for several choices of the stop 
mass. We also show the actual number of t l~ events, taking into account the branching 
ratio for the reaction i t+Na t+l^+l^ +Ni, for an integrated luminosity of £ = 300 fb^^ 
and for a possible upgrade with C = 1 ab^^. For a stop mass below 800 GeV the prospects 
for such a measurement look promising. 



Mi, 




N[ti+£-] 
£ = 300/6-1 Q^-i) 


500 GeV 


300 


7300 (24000) 


800 GeV 


20 


560 (1800) 


1 TeV 


4 


120 (400) 


1.2 TeV 


1 


30 (100) 



Table 1: The production cross-section for iit^ is shown in the middle column. The branching ratio 
for the reaction i ^ t + Na ^ t + 1+ + + Ni was calculated using Mj = 300 GeV, M^^ = 140 GeV, 
Afj^^ = 100 GeV, and assuming that the gluino and squarks are sufficiently heavy to have little effect. 
(Under these assumptions and wino/bino dominated ^2^1 the branching ratio for t tN2 is slightly less 
then 1/3 because of the top's mass, and those for iV2 e^e~Ni or fi^ ^~ Ni are about 1/6 each.) The 
number of t £~ events is then presented in the last column for two different integrated luminosities. 
The effective number of events is doubled if one combines the t£'^£~' and <^£-^+asymmetries. 



In the next section we take into account the experimental difficulties in making such 
a determination. We propose ways of overcoming these limitations and try to evaluate the 
corresponding reduction in signal sensitivity. 
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0.08 



0.06 




0.04 



0.02 



Figure 3: The asymmetry parameter rj is plotted against /ii = A/?^ /-^'^^ ^'^^ several choices of the 
CP-phase A(p. As expected, rj diminishes as the CP-phase decreases. 



5. Experimental limitations 



In this section we consider the degradation of the signal due to experimental limitations. First, 
we address an issue already mentioned above, namely that the triple product is measured in 
the lab frame, and if the LSP escapes detection there is no way to reconstruct the rest frame 
of the stop. In other words, the asymmetry parameter r] was computed in the neutralino 
frame, which cannot be reconstructed. 

Let's imagine an event where the momenta are such that in the neutralino's rest frame 
we have, 

pt • {pi+ xpi-)= p! {pf+p^i- - p\+p1~) > 

This is a contribution to N^. This quantity is still positive even in the stop's rest frame since 
the boost is only along the z-axis and it cannot flip the direction of the top momenta. If the 
stop was produced at rest in the lab frame the signal would be unaltered. However, the stop 
itself is in general boosted with respect to the lab frame. An arbitrary boost can turn this 
contribution to in the stop's frame into a contribution to in the lab frame. It can 
either flip the sign of pt or it can change the transverse orientation of p/+ with respect to pi-. 

The stop - anti-stop pair is produced mainly via gluon fusion and so, owing to the gluon 
distribution function, the stops are produced very close to threshold. However, the overall 
center of mass can be quite boosted with respect to the lab. Therefore, while the stop has very 
little transverse momenta, it does carry a non-negligible momentum along the beam direction. 
In Fig. § we used Pythia |38] to produce a plot of the distribution of stop longitudinal velocity 
in the lab frame for several choices of stop mass. 
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Figure 4: The stop velocity distribution in the lab frame for several choices of stop's mass. All 
distributions are normalized to unit area. As the stop's mass increases the distributions peak at a 
lower /3, in agreement with expectations based on the PDF. 

Therefore, we must account for the possible flip of iV+ into N- (and vice-versa) due to 
the initial boost of the stop. If we denote by w the probability flip then the asymmetry 
parameter rj in the lab frame is given by 



Vlab 



jy(J.ab) _ jy^ilab) 
j^ilab) ^ j^(lab) 



V 



UN) , ^(N) 



(5.1) 



where the dilution factor D is simply 



V = l-2w. 



(5.2) 



One might expect that w ^ as the stop's mass increases since the initial boost is 



remains 



diminished. However, this limit holds true only if the mass difference — m 
fixed. When this mass difference increases, w increases as well. To understand this point, 

increases, all the momenta defining A^4. are on 



Na 



notice that when the difference 
average increased. So, while it is true that it becomes harder to change the sign of pt, it is 
easier to change the orientation of pi+ with respect to pi- . In Fig. |5| we plot the probability 

fixed) as well as a function 



m 



Na 



for a flip, w, as a function of the stop's mass (keeping 
of the mass difference itself (keeping mi fixed). 

From Fig. ^ it is clear that the dilution factor does not present a very serious problem. 
Unless there is a very large splitting between the stop and neutralino masses the probability 
of flip is about w > 0.33. This will translate into a dilution factor of about < 0.1 which 
represents an increase in the number of events needed of about an order of magnitude. 
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Figure 5: In the left pane the probability w of iV+ in the neutralino's rest frame to flip into 7V_ in 



the lab frame (or 7V_ into is plotted as a function of the stop mass keeping 
the right pane, w is plotted as a function of the mass difference, keeping fixed. 



rriff fixed. In 



There are several related issues concerning the efficiency of identifying the top or anti- 
top, determining its charge^ and momentum, determining whether the from the decay is 
associated with the initial t or the (assuming they are pair produced), separating the cascade 
leptons from others in the process, etc. These questions, and how they affect the number of 
identified events, will require an extensive and careful numerical simulation. While we do not 
attempt such a study here but leave it for future research, we would like to make a number 
of comments relevant for such a study. A favorable situation is that in which A'^i and are 
dominantly bino and wino, respectively, and the gluino is heavy. Then the dominant decays 
of the tL should be into b and t N2, with relative rates ~ 2 : 1. When one of the stops 
decays to a chargino and the other cascades, it may optimistically be possible to identify the 
t or and determine its charge by tagging on the charge of the lepton from the chargino 
decay (especially if it is of a different flavor from the dilepton). In this case, one can combine 
the asymmetries from i and i'^ decays; i.e., in analogy to Eq. ( p.9| ) define 



Vs 



N+{i) + N+{i^) + N^{i) + iV_(t^) ' 



(5.3) 



where the 6 angles for the i and i'^ events are defined in the caption to Fig. |^. It is easy 
to show that the theoretical expectation for ijsum is the same as for rj (except for removing 
the strong phase term), while effectively doubling the number of available events, i.e., the 
required (I?^r/^^)~^ is the total number of t l^£^ and t^l^i'^ cascades. 

Events in which the associated top decays semi-leptonically provide another handle on 
its identification and charge, but at the expense of possible additional confusion about which 
lepton is from the top, and missing the neutrino momentum, which renders the reconstruction 



^Averaging over both the top and anti-top processes, which correspond to (T) 
the CP-violating effects and just leaves the strong phase contribution. 



(T) in Section Q, cancels 
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of the top's momentum impossible. None of the observables above really require the top's 
momentum, but only its direction in the lab frame. If the top is highly boosted, it may 
not be necessary to fully reconstruct its momentum to infer its direction in the lab frame. 
More generally, the reduction in signal due to an imprecise determination of the top's direc- 
tion is a detailed numerical question involving an event simulator which we leave for future 
investigation. 

There may also be a non-negligible number of events in which the t decays to t N2 and 
the to t^Ni, or vice versa. Assuming that one determines the charge of t or by its 
leptonic decay, there is still the ambiguity of whether the dilepton is correlated with the 
top or anti-top. One possibility is to simply include the wrong pairing (calculating T (T) 
for the pairing with the t (f^)), and argue that on average it does not contribute to any 
CP violating observable. This effectively sums the asymmetries from the t and t^i~i~^ 
cascades, but without gaining the factor two statistical advantage discussed above in the case 
of one chargino decay. It is possible that the combinatorics can be resolved with a more 
sophisticated analysis involving isolation cuts, energy cuts, etc. For example, if the anti-stop 
decays directly into an anti-top and the LSP, the anti-top is on average more energetic than 
the top coming from the other branch. Ordering the jets according to their pT may reduce 
the combinatorics and help identify the correct pairing. Again, such schemes will require a 
careful numerical simulation. 

Another possibility is to use bottom quarks instead of tops, i.e., consider the reaction 
b ^ b + N2 b + + l~ + Ni. In this case one could reconstruct the full b momentum. 
However, the efficiency of directly determining the b charge is low. These issues could be 
resolved if the opposite side b or b'^ decays to t C~ or C"*", but it may be difficult to 
know whether one began with stops or sbottoms if they are close in mass. For example, 
if one observed t bH'~l'^l~ , the could be associated either with i ^ t i~ ov with 

b'^ — > b'^ i~ . The asymmetries expected from each possibility would be the same for a 
wino-dominated Similar to the discussion of f^Ni above, one could count each event 
twice, once for each possibility, assuming that the wrong pairing does not contribute to the 
asymmetry. The required sum of identified t and b'^ cascades is twice the expression in Eq. 
(|2.12|) , but this number could include the and b cascades if the t%£'~^£~£~^ events are 
combined appropriately. 

A third possibility is to use the asymmetry between quarks and anti-quarks in the par- 
ton distribution function (PDF) of the proton. When considering stop pair-production the 
dominant mode is gluon fusion since there are practically no tops in the proton PDF. In this 
case, stops and anti-stops are produced in equal amounts. However, when considering the 
production of u's and d's there are many more relevant channels. The valence quarks play 
a significant role and associated production {g + q ^ 9 + 0} can dominate. Since valence 
quarks in a proton-proton collider are mostly quarks and not anti-quarks, it is considerably 
more likely to produce a squark than an anti-squark. Therefore, in effect, we know that we 
are observing the reaction and not its CP conjugate. There are several problems with such 
an approach. First, we must include the contribution from all partons (mostly u's and d's). 
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This is in principle easy to incorporate into the present calculation and amounts to a trivial 
addition of the different contributions. The more serious problem is the existence of multijets 
in the event and a reduction of the signal due to combinatorics. 

There is also the problem of the strong phases discussed in Section ^ which in principle 
contribute to the expectation value of the triple-product. One source involves the exchange 
of a photon between the two leptons in the cascade. However, this is of 0(a/7r), and not 
competitive unless the CP-violating phases are small. Another is associated with the phases 
in the slepton propagators from their finite width. This leads to a non-trivial effect in the 
on-shell case for fixed lepton energies, but vanishes when integrated over their energies, as 
commented in the appendix. It is negligible in the more interesting off-shell case. In general, 
the strong phase effect can be eliminated by combining the asymmetries for t and t^, as in Eq. 
( |2.7D or (^]^). However, if one relies on the asymmetry between quarks and anti-quarks in 
the PDF, one cannot form such a combination. In this case, one must rely on the theoretical 
estimate that the strong phases effects are small. 

6. Conclusions 

The triple product and the related asymmetry parameter observable presented in this paper 
are sensitive to CP-phases in the cascade decay of stops, t ^ t + l'^ + l~ + Ni. The phase 
combination that appears in this reaction involves the phase difference between the wino- 
slepton-lepton and bino-slepton-lepton (more generally, A''2 and A^i) couplings. As pointed in 
the introduction, this phase combination is not bounded directly by EDM experiments since 
it does not require a higgsino insertion and is therefore independent of the ^-parameter's 
phase. 

In the case of an on-shell cascade decay, the signal is too small to be observable, requiring 
more than 10^ events to reach experimental sensitivity. However, if the spectrum is such that 
the reaction proceeds via an off-shell slepton, the signal is greatly enhanced. We find that 
about (10^ — 10^)/ sin^(2A(/9) events are needed to constrain the CP-phase Aip. This number 
may improve dramatically if some experimental difficulties discussed in the text are resolved 
or may increase if these turn out to be more severe. At any rate, this number is low enough 
to be taken seriously as a viable observable for probing some of the MSSM's CP-phases at the 
LHC. Rough estimates for the required number of events are given in Table |l[ For large CP 
phases the effect may be observable for stop mass as large as 800 GeV prior to a luminosity 
upgrade and even higher thereafter. A possible luminosity upgrade and a favorable spectrum 
will place the signal well within the experimental sensitivity and help probe a combination of 
the phases which is currently inaccessible via the EDM experiments. 

Our goal has been to illustrate the general possibility and point out the difficulties, not 
to examine the full parameter space. A more systematic study, including a full numerical 
simulation of the events and detector performances, would be very useful. Similar effects 
might also be observable in other channels which may cover different regions of parameter 
space. 



15 



Acknowledgments: We would like to thank G. Kane and T. Han for useful discussions 
in the early stages of this project. The work of L.W. and I.Y. is supported by the National 
Science Foundation under Grant No. 0243680 and the Department of Energy under grant 

# DE-FG02-90ER40542. P.L is supported by the Friends of the IAS and by the NSF grant 
PHY-0503584. The work of G.P. was supported in part by the Department of Energy ^ 
DE-FG02-90ER40542 and by the United States-Israel Bi-national Science Foundation grant 

# 2002272. Any opinions, findings, and conclusions or recommendations expressed in this 
material are those of the author(s) and do not necessarily reflect the views of the National 
Science Foundation. 



A. A derivation of the asymmetry parameter r] in the neutrahno's rest frame 



In this appendix we give the details of the calculation of the asymmetry rj in the neutrahno's 
rest frame (defined in Eq.( |2.9D ). A computation of the expectation value {t^yap Pf^PtP'i+pf- ) 
is a straightforward modification of the derivation below. 

The differential decay width for the reaction is 



2Mi ' ^ ' 



where M is the invariant amplitude. The Feynman rule for the sfermion-fermion-neutralino 
coupling is 




^ {g^ Pl + G{- Pn) . (A-2) 



Our notation follows that of |39|, where the explicit expressions for the G's in terms of the 
mixing matrices can be found. 

The invariant amplitude consists of two parts, corresponding to the Feynman diagrams 
of Fig. H (Th ere are two diagrams contributing to the process as a result of the Majorana 
nature of the neutralinos) . Ignoring the masses of the external leptons and allowing them to 
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have different flavors the two parts are 



ki - Mi 



MPt) {G%*Pl + G'XPn] { -P^] (gI^Pl + Gl-Pn) v{p,,) 



^) ^(PnJ (Gjr^L + G'^f*Pn) v{p,-\ (A-3) 
where q = Pt — Pt, ki = p,+ + p<r , k2 = p,- + Pm , and the masses are taken to be real (i.e., 

5. 1 1 1 



^2 



all the complex phases are absorbed into the couplings). In deriving A-3 we have used some 



of the identities of Appendix C of |33]. 

The asymmetry r] depends on the terms proportional to the Levi-Civita tensor e^^ap. 
Terms in the reaction containing a non-vanishing e^uap can only come from the interference 
terms and must contain 4 independent vectors. We find 



2 Re MaMl D 4 Im (or - a^) e^.^p Pt+ P,- (A-4) 

spin 

1 1 1 1 



X 



_^ Na "a Na ^ Ni, Nf, Nf, 1 /j. 1^ h ^ li k k 

where is given by 

«^ = {GtG%G^GllCGl^^^ (A-5) 
+ M^M^^ {gXG%gXgI-*gI-,G% 



+ c.c 



+ ^Na^N, [g^Ig^-^giicg:]^^^ 

and is simply given by an with L ^ R. In principle, there is also a term proportional to 
Re {an — ai) from the finite width part of the slepton propagators (see Section |^. In the 
fiavor diagonal case this term is proportional to the difference between the Z"*" and l~ energies, 
and vanishes when integrated over phase space. 

If we neglect all the off-diagonal mixing matrix elements, and take the neutralino Na = Nf, 
to be on-shell (g^ = M?. ), the expression simplifies to 
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(A-6) 



where we have defined 

Gt^a'n, Gl-^^gt, = (A-7) 

and similarly for the left handed couplings. Writing each coupling constant as, g = \g\e'^^ 
( |A-6| ) can be expressed in terms of the complex phases i^p as 



Im {aR - ai) 



2M 



Na 



\9r 



(A- 



\9T?) 



bkl'l^k'l'sin 2(^^-^^«) 



The event geometry in the neutralino's rest frame is depicted in Fig. ||. The incoming 
stop and outgoing top define a z-axis with the top pointing in the positive direction. This 
z-axis is in an arbitrary orientation with respect to the lab frame's beam pipe-line axis (since 
the stop is a scalar its decay is isotropic). Momentum conservation in the neutralino's rest 
frame forces the di-lepton and the LSP to lie in the same plane. In other words, the di-lepton 
defines an orthogonal to the plane, 



n=pi+ X pi- 

where h itself is oriented with respect to the z-axis, 

h ■ pt = cos 9 



(A-9) 



(A-10) 



We say that h is in the upper-hemisphere (iV+) if cos0 > or the lower-hemisphere (N-) if 
cos6' < 0. 

A non-zero expectation value for Pf{pi+ xpi- ) translates into a non-zero expectation value 
for — . As far as this difference is concerned the only relevant part of the amplitude is 
the one involving the e^ua/s piece. We are left with evaluating the integral 

f.7r/2 



dPSA 



4Im (or - ol) ef„yaf3 Pi Pt Pi^ Pi 
1 „ / 1 



P 



(A-11) 



Re 



A similar expression holds for only with the limits on the integrals being (tt/2,tt). In the 
neutralino's rest frame Pf = Pt, and therefore 

EiPt • iPi+ xpi-)- Etpi ■ {pi+ X pi_) (A-12) 



U u a 
<^tiuaP Pi Pt Pl+ Pi- 



^N^Pt ■ iPl+ >^Pl~) 

^7V„ \Pt\ \Pi+ \ 1 COS 6* sin ( 
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where cos is defined as above and 4> is the angle between pi+ and pi- . 
The four-body phase space can be written as 

dPS^iPt -^pt+ Pi+ + Pi-- +PnJ = dPS2{pi -^Pt + q) (A-13) 

X dPSz{q^pi+ +pi- +Pff^), 

where the 2-body phase space integral is given by 

The 3-body phase-space integral can be written in terms of dimensionless variables (see for 
example Q) 

Ml 

dPS-i{q^Pi+ +pi- +p^^) = -—^dx+dx- dcos9, (A-15) 

where 

x+ = — — and X- = ——, (A-16) 
and the limits of integration are 

0<a;_<l-/ui, 1 - fii - X- < x+ < 1 - — , (A-17) 

I — X- 

where fii = M^. /M^. . Despite its appearance, the integration domain is symmetric over 
In the narrow-width approximation, the neutralino's propagator is 

^ ^ 6{q'-Ml), (A-18) 



Na' Na Na 

and the dq^ integration can be done trivially, setting = M| everywhere else in the 
expression. Conservation of momentum fixes the angle between pi+ and pi~ to be 

sin</) = 2 {{I — fii — x+ — X- + -|- 2;+ -|- /Ui — 1))"'"''^ • (A-19) 

Noting that = (q—pi-)'^ = (1 — x.) and A;! = iq—Pi+)'^ = (1 — x+), the expression 
for simplifies to 



Q^t a\ I A A 2;+x_sm( 
cos yd cost' / dx+dx- 



cose=o J (1 - a:+ - (l - x_ - ^;-) 



where nj = Mf/M^ 
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Integrating over cos 6, we arrive at the result quoted in the text for the difference N^ — N- 
for a given integrated luminosity C, 



iV+ -N_ 



1 



M 



Na 



2567r3 



dPS-2 

2Mr 



\pt\ 



M 



(A-21) 



Na 



X / dxj^dx- 



(1 - ;Ui - x+ - X-_ + 2;+X-_)(a;+ + x+ + - I] 



1/2 



(1 - X+ - III) (1 - X_ - /X;) 

For the total number of events + we must include the entire amplitude 

iVtotal r dPSi 



2Mr 



\Ma? + \Mb\'' + 2-Re{MaMl)) . 



(A-22) 



Following the same path outlined above it is straightforward to arrive at the result quoted in 
Eqs-dU) and Q. 
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